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^ '. ABSTRACT 

We present the results of a deep (J — 19.1 mag) infrared {ZYJHK) survey over the full 
a Per open cluster extracted from the Data Release 9 of the UKIRT Infrared Deep Sky 
Jl^ ■ Survey Galactic Clusters Survey. We have selected ^ 700 cluster member candidates in 

^ ' ~56 square degrees in a Per by combining photometry in five near-infrared passbands 

^ . and proper motions derived from the multiple epochs provided by the UKIDSS GCS 

DR9. We also provide revised membership for all previously published a Per low-mass 
stars and brown dwarfs recovered in GCS based on the new photometry and astrometry 
provided by DR9. We find no evidence of if-band variability in members of a Per with 
dispersion less than 0.06-0.09 mag. We employed two independent but complementary 
\ methods to derive the cluster luminosity and mass functions: a probabilistic analysis 

^\ • and a more standard approach consisting of stricter astrometric and photometric cuts. 

\^ ' We find that the resulting luminosity and mass functions obtained from both methods 

are consistent. We find that the shape of the a Per mass function is similar to that 
^ ■ of the Pleiades although the characteristic mass may be higher after including higher 

I mass data from earlier studies (the dispersion is comparable). We conclude that the 

mass functions of a Per, the Pleiades, and Praesepe are best reproduced by a log- 
normal representation similar to the system field mass function although with some 
variation in the characteristic mass and dispersion values. 

, Key words: Techniques: photometric — stars: low-mass, brown dwarfs; stars: lu- 

■ minosity function, mass function — galaxy: open clusters and associations: individual 

(Alpha Per) — infrared: stars 
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1 INTRODUCTION 

The shape of the Initial Mass function (IMF) is of prime 
importance to understand the processes responsible for the 
formation of stars and brown dwarfs. The d efinition and th e 
first estimate of the IMF was presented in ISalpete 3 (Il955ll . 
Our knowledge of the IMF has now improved both at the 
high-mass and low-mass ends. The mass spectrum in open 
clusters and in the field, defined as dN/dMocM"" {a is the 
exponent of the power law and equivalent to 1, where x is 
the slope of the logarithmic mass function), is currently best 

* Based on observations made with the United Kingdom Infrared 
Telescope, operated by the Joint Astronomy Centre on behalf of 
the U.K. Particle Physics and Astronomy Research Council, 
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fit by a three segment power law with a = 2.7 for stars more 
massive than 1 M©, a = 2.2 between 1 and 0.5 M^, and a 
= 1.3±0.5 in the 0.5-0.08 Mq mass range (iKroupal |2003 ) . 
Alternatively, a log-normal function with a cha racteristic 
mass around 0.2-0.25 M© and dispersion ~ 0.55 l|Chabrieii 
[2003', '2005') provides a good match to current observations 
for the system mass function in the field. The advent of large- 
scale optical and near-infrared surveys towards open clusters 
extended the mass spectrum to the substellar regime but a 
consensus has yet to emerge on the detailed shape. 

a Per is one of the few open star clusters within 200 
pc of the Sun and younger than 200 Myr. The cluster 
is located to the north-east of the F5 V supergiant Alpha 
Perse i at a distance o f ^17 5-190 pc (jPinsonneault et al.l 
1 19981 : iRobichon etHI 1 19991 ) with a revised distance of 
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172.4±2.7 pc from the re-reduction of the Hipparcos data 
llvan Leeuwerj|2009l). The cluster members have solar metal- 
licity (|Boesgaard fc Friell[l990l ) and the extinction along the 
line of sight is estimated as Ay = 0.30 mag with a possi- 
ble differential extinction (|Prossei|[l993 ). It has been well 
studied, though less frequently than the Pleiades due to 
a smalle r proper motion (( /Xg cos (5,^t3) = (-1-22.73,-26.51) 
mas/yr; Ivan Leeuwen|[200g ) and a much lower galactic lat- 
itude (b = —7° vs. —24°). Despite being further away than 
the Pleiades (170 pc vs. 120 pc), a Per is a good target for 
substellar studies because it is younger than the Pleiades 
(85±10 Myr vs. 125±8 Myr), placing the lithium depletion 
boundary at 7c ~ 17.7-17.8 mag for both clusters. 

Multi-wavelength surveys and spectroscopic follow-up 
observations have been performed in a Per to extract a clean 
sequence of cluster members from high-mass stars down to 
brown dwarfs. The fir st proper motion surve y in the clus- 
ter was performed bv lHeckmann et al.l |l95^ and comple- 
mented by photometry from lMitchelll (|l960l ). yielding about 
60 probable me mbers (HE obje cts) whose final membership 
was revised by IProsseiJ l| 19921 ). The me mbership of addi- 



tional candidates pro posed bv iFresneaul ([l98fl) was subse- 
quently established bv lProsseil (Il992h. Lower mass members 
(AP sources) were e xtracted bv IStauffer et all (|l985l ) and 
IStauffer et aD (|l989l ) on the basis of th eir prop e r mot ion, 
photometry, and spectral characteristics. iProsseJ (| 19921 ) ex- 
amined the Palomar photographic plates to extract new low- 
mass proper motion and photometric members down to a 
spectral type of M4 over a 6° by 6° field. Additional low- 
mass photometric candidates wer e reported fro m a deeper 
optical survey in a smaller area (|ProsseJll994l') as well as 
from X-rays observ a tions with ROSAT ( Randich et al.lll9^: 
Prosser et al ] 1 19961 : IProsser fc Randichlll998l : IProsser et all 



19981 ). The first brown d warf candid a tes w ere spectro- 



scopically confirmed by iStauffer et al.l (| 19991 ). yielding a 
lithium age of 90±10 Myr, twice the turn-off main-sequence 
age (50 Myr; iMermilliodI Il98ll ). A revised value of the 
age derived from the lith i um m ethod was published by 
iBarrado v Navascues et al.1 l|2004l ). estimated to 85±10 Myr. 
A deep optical survey complemented by near-infrared pho- 
tometry extended the cluster se quence down to 0.03 M0 
ijBarrado v Navascues et al.ll200j ). The best fit of the slope 
of the mass function was obtained for a power law index a 
= 0.59±0.05 over the 0.3-0.035 M© m ass range, in agree- 
ment with estimat es in the Pleiades l|Dobbie et al.l I2OO2I : 
iMoraux et ai]|2003l ) at that time. A wider survey based on 
photographic plates bv lDeacon fc Hamblvl I2OOJ) derived a 
power law index a of 0.86 (0.67-1.00) over the 1.0-0.2 M© 
range from a sample of high probabil i ty me mbers over ~250 
square degrees. Finallv. [Lodieu et al.l (|2005l ) extracted about 
20 new infrared photometric candidates from a deep AT-band 
survey of 0.7 square degree prev iously covered in the opti- 
cal bv iBarrado v Navascues et al . (2002 ). Additionally, 24 
probable candidates from lBarrado v Navascues et al.l ()2002l ) 
were confirmed as spectroscopic members with masses be- 
tween 0.4 and 0.12 M©. 

The UKIRT In frared Deep Sky Survey (UKIDSS; 
iLawrence et l2007^ is a deep large-scale infrared sur- 
vey conducted w ith the wide-field camera WFCAM 
ijCasah et al.ll2007l ) on UKIRT (Mauna Kea, Hawai'i). The 
survey is subdivided into 5 components: the Large Area 
Survey, the Galactic Clusters Survey (hereafter GCS), the 



Galactic Plane Survey, the Deep Extragalactic Survey, and 
the Ultra-Deep Survey. The GCS aims at covering ~1000 
square degrees in 10 star-forming regions and open clusters 
down to K = 18.4 mag at two epochs. The main scientific 
driver of the survey is to study the IMF and its dependence 
with environment in the substellar regime using an homo- 
geneous set of low-mass stars and brown dwarfs over a large 
area in several regions. 

In this paper we present the a Per mass function over 
~56 square degrees derived from the UKIDSS GCS Data 
Release 9 (DR9). This is the second paper of its kind after 
the an alysis of the Pleiades cluster presented in lLodieu et al.l 
(i2012h . In Section [2] we present the photometric and astro- 
metric dataset employed to extract member candidates in a 
Per. In Section |3] we review the list of previously published 
members recovered by the UKIDSS GCS DR9 and revise 
their membership. In Section |4] we outline two methods for 
deriving the cluster luminosity function. One method relies 
on a relatively conservative photometric selection followed 
by the calculation of formal membership probabilities based 
on object positions in the proper motion vector point dia- 
gram (Section l4.1|l . The second method applies a more strin- 
gent colour cut followed by an astrometric selection based on 
the formal errors on the proper motions for each photomet- 
ric candidate compared to that of the cluster (Section I4.2|l 
for which we test the level of contamination (Sect.[5)|. In Sec- 
tion |6] we discuss the K-hasiA variability of cluster member 
candidates in a Per. In Section[7]we derive the cluster lumi- 
nosity and (system) mass function and compare it to other 
clusters studied as part of the GCS (Pleiades and Praesepe), 
and the field. 



2 THE UKIDSS GCS IN a Per 

The UKIDSS GCS DR9 released ^56 square degrees ob- 
served in five passbands {ZY JHK; iHewett et aLlIiooil ) in 
the a Per open cluster over a region defined by RA=44-60 
degrees and dec=44-54 degrees. 

We have selected all good quality point sources in a Per 
detected in at least JHKl (where Kl stands for the first 
Tf-band epoch) and, where available, in Z, Y, and K2 (sec- 
ond A"-band epoch). We imposed a request on point sources 
only in JHK and pushed the completeness towards the faint 
end by imposing limits on the ClassStat parameters (be- 
tween —3 and -f3) which classify the point-likeness of an 
image. The Structured Query Language (SQL) query used 
to select sources along the line of sight of the a Per is iden - 
tical to the query used for the Pleiades ijLodieu et al.ll201^ ). 
The SQL query i ncludes the cross-ma tches with 2MASS 
(jCutri et al.ll2003l : ISkrutskie et aLllioO^ to compute proper 
motions for all sources brighter than the 2MASS 5a com- 
pleteness limit at J = 15.8 mag as well as the selection of 
proper motions from multiple epochs provided by the GCS. 
We used the GCS proper motion measurements in this work 
as they are more accurate due to the homogeneous cover- 
age, completeness, and spatial resolution of the UKIDSS im- 
ages and the detailed rela tive astrometric mapping employed 
(jCollins fc Hamblvl20l3 ), and of course the GCS proper mo- 
tions are available for objects that are too faint for 2MASS. 
We limited our selection to sources fainter than Z = 11.6, 
Y = 11.4, J = 11.0, H = 11.5, Kl = 10.0, K2 = 10.4 mag 
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Figure 1. Completeness of the GCS DR9 dataset in the a Per 
cluster in each of the six filters. The polynomial fit of order 2 is 
shown as a red line and defines the 100% completeness limit of 
the GCS DR9 in each passband. 

to avoid saturated point sources. The completeness limits, 
taken as the magnitude where the straight line fitting the 
shape of the number of sources as a function of magnitude 
falls off, are Z = 20.0, Y = 19.6, J = 19.1, H = 18.4, Kl = 
17.6, and K2 = 18.1 mag (Fig.©. 

The query returned 2,643,045 sources with J — 11.0- 
21.2 mag over ~56 square degrees towards the a Per clus- 
ter. The full coverage is displayed in Fig. [2]and the resulting 
{Z — J,Z) colour-magnitude diagram is shown in Fig.Oalong 
with previously published member candidates (black filled 
dots). Note that theoretical isochrones plotted in this paper 
were specifically computed for the WFCAM set of filters at 
an age of 90 Myr (downloaded from France AUard's web- 
page jj. We combined the NextGen and DUSTY isochrones 
for effective temperatures above and below 2700 K, respec- 
tively, to convert magnitudes into masses (Section [7|). 



3 CROSS-MATCH WITH PREVIOUS 
SURVEYS 



There are 455 probable members known in a Per ex- 
tracte d from previous prope r motion and o p tical sur- 



veys ( 


Heckmann et al.1 Il95d;_ Mitchelll 1960l: Fresneaul 


198d: 


Stauffcr ct al. 19851. \l98^: Prosseij 19921. 19941: 


Prosser & Randichil99S: iProsser et al. 


19981: Stauffer et al.1 


19991: Barrado v Navascues et al.l 20021; 


Lodieu et al. 2005h. 



and an addition al 300 high-probabi lity (p J5 60%) member 
candidates from iDeacon fc Hamblvl (|2004 V 



We cross-matched catalogues from earlier studies with 
our full sample of over ~2.5 million sources retrieved from 

^ Prance AUard's Phoenix web simulator can be found at 
|http://phoenix.ens-lyon.fr/simulator/index.faces| 



APer GCS DR3 coverage 




J , , , \ I i , I , i i_ 

4 2 

RA offset (degrees) 



Figure 2. Coverage from the UKIDSS GCS DR9 in the a Per 
open cluster in the standard angular plane coordinates (^, rj) 
choosing (ra,dec) = (51, 49) degrees as the cluster centre. The 
total area covered is about 56 square degrees. The holes present 
in the coverage are due to the rejection of some tiles after quality 
control. GCS DR9 member candidates identified in this work are 
overplotted as black filled dots. 
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Figure 3. {Z — J,Z) CMD for ~56 square degrees in the a Per ex- 
tracted from the UKIDSS Galactic Cluster Survey Data Release 
9. Previously published member candidates in a Per are overplot- 
ted as filled dots. The mass scale is shown on the right hand side 
of the diagrams and extends down to 0.03 Mq, according to the 
NextGen and DUST Y models assuming an age of 90 Myr an d a 
distance of 172.4 pc l|Baraffe et al.lll998l : IChabrier et~alll2000ll . 
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Table 1. Updated membership of member candidates identified in a Per by ea rlier studies and recovered in the GCS DR9 sample. Papers 
study i ng a P er over the pa s t decades and considered in this work are: iH eckmann ct al. 1956; Mitchell 1960; Frcsncau 1980; Stauffcr ct al] 
1985, 1989; Prosscr 1992, 1994; Prosscr fc RandichI [1991: iProsser et al.lll998l: IStauffer et al.i Il999l : IBarrado v Navascues et al.l I200a 
Deac on & Hambly 2004; Lodieu ct al. 2005). Columns 2 and 3 give the numbers of sources published by the reference given in Col- 
umn 1 and the numbers of sources recovered in GCS DR9, respectively. Column 4 (named No_DR9) is subdivided into several columns 
to give the reasons why some of the sources from earlier studies are not covered: "Bright" stands for objects brighter than the GCS 
saturation limits, "Outside" stands for sources outside the GCS DR9 coverage, "No_mag" stands for sources missing at least one of the 
J, H, or K magnitudes, ">3"" stands for sources beyond the 3 arcsec matching radius used in our study, and "Flag" stands for sources 
whose GCS flags arc too bad to be included in our catalogue of point sources. Columns 5 and 6 give the numbers of high-probability 
members (p^^40%) and non members (NM) according to our probabilistic approach (first number) and method #2 (second number). 
The last column gives the percentages of sources recovered in the GCS DR9 (ratio DR9/A11). 
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GCS DR9 to locate the cluster sequence in various colour- 
magnitude diagrams. We recovered a total of 426 known 
members in a Per after removing multiple detections present 
in various catalogues (Table IXTj) . The numbers and percent- 
ages in brackets in the second and sixth column of Table [T] 
consider previously published sources lying in the magnitude 
range probed by the GCS. We also made a detailed analysis 
of the 629 previously known members not recovered by our 
SQL query. The numbers are given in the fourth column of 
Table lAll which is divided into five sub-columns. Most of 
these sources are either missing an image in J, H, or Kl or 
are not covered by the GCS (223 or 35.5%) or are brighter 
than the saturation limits set in our query (205 or 32.6%) 
or are very likely proper motion non members (48 or 7.6%). 

4 NEW SUBSTELLAR MEMBERS IN a Per 
4.1 Probabilistic approach 

4.1.1 Method 

In this section we outline the probabilistic approach we em- 
ployed to select low-mass stars and brown dwarf member 
candidates in a Per using photometry and astrometry from 
th e UKIDSS GCS DR9. Th is me thod is described in detail 
in iDeacon fc Hamblvl (|2004l ) and iLodieu et all |20o3). The 
main steps are: 

(i) Define the cluster sequence using candidates published 
in the literature within the area covered by the GCS DR9 

(ii) Make a conservative cut in the [Z — J ,Z) diagram 
to include known members and any new cluster member 
candidates defined as (.^ ^ 16.5 and Z ^(11.5-1- 5. Ox (Z-J)) 
OR (Z ^ 16.5 and Z ^(8.5 + 8.0 x (Z-J)) displayed as solid 
black lines on the top-left panel in Fig. [S] 



(iii) Analyse the vector point diagram in a probabilistic 
manner to assign a membership probability for each photo- 
metric candidate with a proper motion measurement (Sec- 
tion |4T2l). 

(iv) Derive the luminosity and mass function by summa- 
tion of membership probabilities to provide a statistically 
complete sample. 



4.1.2 Membership probabilities 

In order to calculate formal memb ership probabilities we 
have used the same technique as in iLodieu et all (|2012l ) to 
fit distri bution functions to proper motion vector point di- 
agrams (|Hamblv et al.|[l995l ). The techn ique differs slightly 
from the o riginal method pres ented in iDeacon fc Hamblvl 
(2004) and ILodieu et al.l (|2007l ) in that the value of sigma 
(a) is fixed by the formal astrometric errors propagated from 
the centroiding errors given by the source extraction soft- 
ware employed upst ream of the WFCAM Science Archive 
l|Hamblv et al.ll2008l ). the main repository of UKIDSS data. 

We refer the reader to the above cited papersfor more 
details and additional equations. First, we rotated the vec- 
tor point diagram so the cluster lies on the y-axis, assum- 
ing a proper motio n of (22.73,-26.51) mas/yr for a Per 
(|van Leeuwen|[2009l ) after applying a very conservative pho- 
tometric selection in the (Z — J,Z) colour-magnitude dia- 
gram. We also note that we used the following rotation for 
the vector point diagram: 

• fj.^' =cos(0.77xPI) X - sin(0.77xPI) x fiy 

• Illy = cos(0.77xPI) X + sin(0.77xPI) x iiy 

We have assumed that there are two contributions 
to the total distribution (fi{fj,x , ^ly) , one from the cluster, 
4>c{lix , liy) , and one from the field stars, (f)f{^x,l-'-v)- The 
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Table 2. Summary of the results after running the programme to 
derive membership probabiUties. For each Z magnitude range, we 
list the number of stars used in the fit (Nb), the field star fraction 
f, and parameters describing the cluster and field star distribution. 
Units are in mas/yr except for the number of stars and the field 
star fraction f. The cluster star distribution is described by the 
mean proper motions in the x and y directions (/i^^ and Hy^) 
and a standard deviation a. Similarly, the field star distribution 
is characterised by a scale length for the y axis (t), a standard 
deviation T,x, and a mean proper motion in the x direction {/ixf)- 
Note that the value of sigma (ct) is fixed by the formal astrometric 
errors. 
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Figure 4. Vector point diagram showing the proper motion in 
right ascension (x-axis) and declination (y-axis) for previously 
known member candidates recovered by the GCS DR,9 (black 
dots) and the new member candidates selected with method #2 
(red dots). 

fitting region was delineated by —50 < jj,x < 50 mas/yr and 
— 50 < /ly < 50 mas/yr. These were added by means of a 
field star fraction /. 

We characterised the cluster distribution as a bivari- 
ate gaussian with a single standard deviation a and mean 
proper motion values in each axis and Hyc- The field 
star distribution was fitted by a univariate gaussian in the x 
axis (with standard deviation T,x and mean fJ,xf) and a de- 
chning exponential in the y axis with a scale length r. The 
use of a declining exponential is a standard method (e.g. 
Ijones fc^Stauffer 1991jl and is justified in that the field star 
distribution is not simply a circularly-symmetric error dis- 
tribution (i.e. capable of being modelled as a 2d Gaussian) - 
rather there is a prefered direction of real field star motions 
resulting in a characteristic velocity ellipsoidal signature, i.e. 
a non-Gaussian tail, in the vector point diagram. This is best 
modelled (away from the central error-dominated distribu- 
tion) as an exponential in the direction of the antapex (of 
the solar motion). 

The best fitting set of para meters were chosen using a 
maximum likelihood method (^see lDeacon fc Hamblvll2004 ). 
However in a deviation from this method we did not fit for 
the standard deviation of the cluster proper motions (a). 
Instead we calculated the mean astrometric error for all ob- 
jects in each magnitude range and used this as our clus- 
ter standard devia t ion. T his fitting process was tested by 
iDeacon fc Hamblvl (|2004l ) where simulated data sets were 
created and run through the fitting process to recover the 
input parameters. These tests produced no significant off- 
set sjn_the_2axametovaJ^ (see Table 3 and Appendix A 
of|D eacon fc Hamblvl [2003 . for results and more details on 
the procedure). Hence, we calculated the formal membership 
probabilities as. 



+ (1 - /)0c ' ' 

We split the sample into eight intervals of magnitudes 
because astrometric errors are a function of magnitude and 
also to improve the contrast between the field stars and the 
cluster. Each band is one magnitude wide and was fitted 
wi th all seven parameters in the same way as described 
in iDeacon fc Hamblvl (|2004l ). There was no fit possible for 
the 20-21 magnitude bins because of the small number of 
sources in this bin. A summary of the fitted parameters from 
the probabilistic analysis described above is given in Table 


4.1.3 Probabilistic sample 

The probabilistic approach yielded a total sample of 10,176 
sources with membership probabilities assigned to each of 
them. This sample contains 728 sources with membership 
probabilities higher than 40% (including known ones previ- 
ously published) listed in Table IBll Tightening this proba- 
bility threshold to 50% and 60% yields samples of 573 (~27% 
less) and 431 (~69% less) member candidates in a Per, re- 
spectively. These high-probability members are displayed in 
Fig. [S] with previously published candidates in a Per plotted 
in black. 

4.2 Photometry and proper motion selection 

In this section we outline a more widely used method (ref- 
ered to as method #2 in the rest of the paper) that we ap- 
plied to select low-mass and substellar member candidates in 
a Per. This procedure consist of selecting cluster candidates 
by applying proper motion selection followed by strict pho- 
tometric cuts in various colour-magnitude diagrams. This 
alternative method provides and independent test of the 
probabilistic approach presented in the previous section. 

The first step was to select all sources with formal er- 
rors on the proper motion within 3a of the mean proper mo- 
tion of the cluster (Fig. |4]), yielding a completeness better 
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Figure 5. Colour-magnitude diagrams showing the member candidates previously reported in a Per (black dots) and all candidates 
extracted from our probabilistic analysis, including known ones (red dots). Upper left: (Z — J,Z); Upper r ight: {Z — K,Z); Lower left: 
(Y - J,Y); Lower r ight: (J - K,J). Overplotted are the 90 Myr NextGen (solid line: lBaraffe et al.lll998h and DUSTY (dashed hne; 
IChabrier et al.ll200(]l l isochrones shifted to a distance of 120 pc. The mass scale is shown on the right hand side of the diagrams and 
spans 0.60—0.03 Mq, according to the 90 Myr isochrone models. The solid black lines in the upper left diagram represent our conservative 
photometric cuts used for the probabilistic approach. 



than 99% assuming normally distributed errors. The main 
advantage of this method is that it does not rely on a sin- 
gle radius for the proper motion selection but rather takes 
into account the increasing uncertainty on the proper mo- 
tion measurements between the GCS epochs with decreasing 
brightness. 

Secondly, we plotted several colour-magnitude diagrams 
(Fig. [5} to define a series of lines based on the position of 
known a Per members identified in earlier studies and pub- 



lished over the past decades (Table [T|. Those lines detailed 
below are plotted in Fig. |6] and improve on the pure proper 
motion selection. We note t hat those criteria a re similar to 
those used for the Pleiades (|Lodieu et al.|[20l3 ) because the 
younger age of a Per compared to the Pleiades is compen- 
sated by its larger distance. 

• {Z - J,Z) = (0.60,12.0) to (1.20,16.5) 

• {Z - J,Z) = (1.20,16.5) to (2.00,20.0) 
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Figure 6. Same as Fig. [5] but only for member candidates in a Per selected using method #2. The YJHK and JHK-only detections 
have been added too for completeness. 



• 




K,Z) = 


(1.20,11.5) to (1.95,17.0) 


• 


(z- 


K,Z) = 


(1.95,17.0) to (4.00,21.5) 


• 


(Y- 


JX) = 


(0.30,11.5) to (0.55,16.5) 


• 


(Y- 


JX) = 


(0.55,16.0) to (1.40,20.5) 


• 




K,J) = 


(0.75,11.0) to (0.75,16.5) 


• 


(J- 


K,J) = 


(0.75,16.5) to (1.70,19.0) 



This selection returned a total of 685 low-mass stars and 
brown dwarfs with Z magnitude ranging from 12 to 21.5, in- 
cluding known ones recovered by the GCS (Table [Cl). This 
total number is similar to the number of high probability 
member candidates — 728 (431) with p > 40 (60)% — in a 
Per identified via the probabilistic approach. 



4.3 Search for lower mass members 

In this section we search for fainter and cooler substellar 
members in a Per by dropping the constraint on the Z-band 
detection and later the Z + Y bands. 

4.3.1 YJHK detections 

To extend the a Per cluster sequence to fainter brown dwarfs 
and cooler temperatures, we searched for potential candidate 
members undetected in Z. We imposed similar photometric 
and astrometric criteria as those detailed in Section [4.21 but 
analysed Z drop-outs as follows: 
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• y ^ 18 and J s; 19.1 mag 

• Candidates should lie above the line defined by {Y — 
J,Y) = (0.55,16.0) and (1.40,20.5) 

• Candidates should lie above the line defined by (J — 
K,J) = (0.75,16.5) and (1.70,19.0) 

• The position on the proper motion vector point diagram 
of each candidate should not deviate from the assumed clus- 
ter proper motion by more than 3cr 

This selection returned 13 additional member candidates in 
a Per (Table lPTj) . All but four of them are indeed undetected 
in the Z-band images and look well detected in the other 
bands after checking the CCS DR9 images. Thus we are left 
with nine bona-fide member candidates. 

4.3.2 JHK detections 

We repeated the procedure described above looking for Z 
and Y non detections. We additionally applied the following 
criteria: 

• J = 18-19.1 mag 

• Candidates should lie above the line defined by (J — 
K,J) = (0.75,16.5) and (1.70,19.0) 

• The position on the proper motion vector point diagram 
of each candidate should not deviate from the assumed clus- 
ter proper motion by more than Scr 

This query returned 36 new candidate members in a 
Per. After checking the CCS images, we retained only eight 
of them as bona-fide candidates because the others actu- 
ally appear in the Z and/or Y images (although detections 
are not reported in the CCS DR9 catalogue) or have no Z 
and/or Y images. The reasons for the rejection of 28 of the 
36 candidates is given in the last column of Table IDll 



5 ESTIMATION OF THE CONTAMINATION 

In this section we estimate the level contamination present 
in our photometric and astrometric selection (method #2). 

The number density of field objects in our final list of 
candidates as a function of mass is obtained in a similar way 
as in Boudreault et al. (2012, submitted to MNRAS). We 
obtained the radial profile of our cluster candidates in three 
mass ranges: above 0.3 M0, between 0.072 and 0.3 Mq, and 
below the hydrogen-burning limit at 0.072 Mq (Fig. |8}. 

However, considering the incomplete coverage of the 
UKIDSS GCS DR9 towards a Per (holes present in the cov- 
erage due to quality control, see Fig. [Sjl , all datapoints must 
be considered as lower limits: we are only par tly complete up 
to the tidal radius of a Per at 2.91° (9.7 pc; iMakarovl [20061 ') 
and up to 3.5° (95% complete in coverage). Consequently, 
the estimated contamination represent an upper limit. 

We used only the number of objects between 3 and 3.5° 
(outside the estimated tidal radius) at each mass range to 
obtain an upper limit of contamination. This gives 2.92 ob- 
jects per square degree for candidates with masses above 
0.3 Mq, 1.57 between 0.072 and 0.3 Mq, and 0.62 objects 
per square degree for our substellar candidates. Within 3° 
from the cluster centre, this gives a contamination of 35.1%, 
15.9% and 50.6% for the same mass range respectively, or 
26.3% for the whole a Per sample within 3°. This level of 
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Figure 7. Difference in the K magnitude {K1-K2) as a func- 
tion of the Kl magnitude for all member candidates in a Per 
selected with method #2. The YJHK and J H K-only detections 
have been added too (dots with open squares and open trian- 
gles, respectively). Typical error bars on the KI—K2 colours as a 
function of magnitude arc displayed as dotted lines. 

contamination brings into agreement within a factor 
of two the luminosity functions derived from both 
selection methods highlighted in this paper (left- 
hand sidepanel of Fig. 

These numbers appear quite large. We stress again that 
these are upper limits, since the coverage is not complete. 
However, we can claim an completeness higher than 90% for 
our cluster candidate list and the determination of our mass 
function. This is justified by the fact that our astrometric 
selection includes all objects within 3a of the cluster's mean 
proper motion (completeness of >99%) and that the lines 
used in our photometric selection go at least 2a bluer from 
the cluster main sequence in all the colour-magnitude dia- 
grams used for the photometric selection (completeness of 
-95.4%). 

Most of the contaminants of our cluster candidates with 
masses above 0.1 Mq would be Galactic disk late-type and 
giant stars, while most of the contaminants of candidates 
less massive than 0.1 Mq would include Galactic disk late- 
type and giant stars, but also unresolved galaxies. 



6 VARIABILITY AT 90 MYR 

In this section we discuss the isT-band variability of the low- 
mass stars and brown dwarfs in a Per using the two epochs 
provided by the GCS. First we considered the candidates 
extracted with method #2, several of them being already 
published in the literature (Tables TAIII . 

Figure [7] shows {K1-K2) versus Kl for all candidate 
members in a Per from method #2. The brightening in the 
Kl = 11-12 mag range is due to the difference in depth be- 



The Alpha Per astrometric and photometric mass function 9 




1 2 3 4 5 6 7 
Distance from cluster centre [Deg] 




1 2 3 4 5 6 7 
Distance from cluster centre [Deg] 



T" 



M < 0.072 



1 2 3 4 5 6 7 
Distance from cluster centre [Deg] 



Figure 8. Radial density plots of our candidate members of a Per 
in three mass ranges: above 0.3 Mq (top panel), between 0.72 
and 0.3 Mq (middle panel), and below the stellar/substellar limit 
at 0.072 Mq (low panel). The error bars on each datapoint are 
Poissonian arising from the number of objects in each bin. The 
dotted horizontal line is the estimated contamination per square 
degree for each mass range. 



tween the first and second epoch, around 0.5 mag both in 
the saturation and completeness limit. This is understand- 
able because the exposure times have been doubled for the 
second epoch with relaxed constraints on the seeing require- 
ment and weather conditions. We excluded those objects 
from our variability study. Overall, the sequence is very well- 
defined and very few objects appear variable in the K-h&nd. 

We selected variable objects by looking at the standard 
deviation, robustly estimated as 1.48 x the median absolute 



deviation which is the median of the sorted set of absolute 
values of deviation from the central value of the K1-K2 
colour. We identified one potential variable object in the 
Kl — 11-12 and 12-13 mag range with differences in the 
A'-band larger than 3cr above the standard deviation. No 
additional variable source was picked up beyond 3a down 
to Kl — 16.5 mag. The candidate selected in the brightest 
bin appears saturated in the second epoch image, suggesting 
that the variability may be caused by the inaccurate pho- 
tometry derived from saturated sources. The other source 
does not look saturated: its variability may be attributed 
to the presence of a faint companion located south-east at 
~1.2 arcsec best visible in the K2 image due to the greater 
depth of the second epoch. This variability analysis is not 
feasible for Kl ^ 16.5 mag due to the small number of a 
Per candidate members beyond that magnitude range. 

We conclude that the level of variability at 90 Myr 
is small, with standard deviations in the 0.06-0.09 mag 
range, suggesting that it cannot account for the dispersion 
in the cluster sequence. The same conclusions are drawn 
from the high probability samp le and are consistent with 
our analysis of the Pleiades (Lodieu et al.| [2012) and Prae- 
sepe (Boudreault et al. subm) samples although we should 
point out that a handful of members are found to be vari- 
able. 



7 THE LUMINOSITY AND INITIAL MASS 
FUNCTIONS 

In this section we discuss the cluster luminosity and mass 
functions derived from the samples of member candidates in 
a Per extracted from both methods described in the previ- 
ous section. We did not attempt to correct the mass function 
for binaries, hence, we compare our results to "system" mass 
functions. Note that contrary to our work in the Pleiades 
and Praesepe, we are unable to estimate the substellar mul- 
tiplicity due to larger scatter in the single-star and binary 
sequences due to crowding. 



7.1 Age and distance of a Per 

Age determinations i n open clusters can vary by up 
to a factor of two (|jeffries fc Navloil 1200 ll ): fitting of 
the upper main-sequence and giant branch (|Mermilliodl 
Il98lf ) compari ng with models including some convec- 



tive overshoot (|Maeder fc MermilliodI 
younger ages than the lithium test 
In the case of a Per, the former 



T983)_tend_to^ld 



( Rebolo et~ai]|l992l ). 

method gives 51 
Myr whereas the latter suggests an a ge between 85±10 
Mvr (iBarrado v Navascues et al.l |2004| '1 and 90±10 Myr 
(|Stauffer et al.l 19991 ). A similar discrepancy has been ob- 
served for the Pleiades ( 77 Myr vs. 120 Myr; Mer milliodI 
Il98ll : IStauffer et al]ll998l ). Moreover, iMevnet et al.l (jl993l ) 
revised the ages of 30 galactic open clusters based on an up- 
dated set of solar-metallicity isochrones (at that time) taking 
into account mass loss and moderate overshooting, yielding 
52 Myr and 100 Myr for the a Per and the Pleiades, re- 
spectively. The latter age for the Pleiades is favoured by the 
fit ting technique of the main-sequence evolution developed 
bv lNavlorl (|2009l ) which quoted a 68% confidence interval of 
104-117 Myr (mean value of 115 Myr), in agreement with 
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Table 3. Values for the luminosity and mass functions (both in linear and logarithmic scales) per magnitude and mass bin for the a 
Per open cluster from the probabilistic approach. We assumed a distance of 172.4 pc and employed the NextGen and DUSTY 90 Myr 
theoretical isochrones. 



Mag 


range 


Mass 


range 


Mid-mass 


dN 


errH 


errL 


dN/dM 


errH 


errL 


dN/dlogM 


errH 


errL 


12.0- 


-12.5 


0.7380- 


-0.6420 


0.6900 


6.01 


3.60 


2.40 


62.60 


37.50 


25.00 


2 


00 


0.47 


0.51 


12.5- 


-13.0 


0.6420- 


-0.5750 


0.6085 


19.49 


5.50 


4.39 


290.90 


82.07 


65.47 


2 


61 


0.25 


0.25 


13.0- 


-13.5 


0.5750- 


-0.5070 


0.5410 


48.33 


8.01 


6.93 


710.74 


117.73 


101.97 


2 


95 


0.15 


0.15 


13.5- 


-14.0 


0.5070- 


-0.4200 


0.4635 


63.98 


9.05 


7.98 


735.40 


103.97 


91.76 


2 


89 


0.13 


0.13 


14.0- 


-14.5 


0.4200- 


-0.3260 


0.3730 


66.16 


9.18 


8.12 


703.83 


97.66 


86.37 


2 


78 


0.13 


0.13 


14.5- 


-15.0 


0.3260- 


-0.2440 


0.2850 


89.67 


10.51 


9.46 


1093.54 


128.16 


115.32 


2 


85 


0.11 


0.11 


15.0- 


-15.5 


0.2440- 


-0.1830 


0.2135 


77.31 


9.84 


8.78 


1267.38 


161.23 


143.91 


2 


79 


0.12 


0.12 


15.5- 


-16.0 


0.1830- 


-0.1390 


0.1610 


70.11 


9.42 


8.36 


1593.41 


214.04 


189.96 


2 


77 


0.13 


0.13 


16.0- 


-16.5 


0.1390- 


-0.1085 


0.1237 


51.75 


8.25 


7.18 


1696.72 


270.35 


235.29 


2 


68 


0.15 


0.15 


16.5- 


-17.0 


0.1085- 


-0.0869 


0.0977 


50.93 


8.19 


7.12 


2357.87 


379.11 


329.58 


2 


72 


0.15 


0.15 


17.0- 


-17.5 


0.0869- 


-0.0703 


0.0786 


19.14 


5.46 


4.35 


1153.01 


328.90 


261.82 


2 


32 


0.25 


0.26 


17.5- 


-18.0 


0.0703- 


-0.0591 


0.0647 


8.26 


4.00 


2.83 


737.50 


357.29 


252.70 


2 


04 


0.40 


0.42 


18.0- 


-18.5 


0.0591- 


-0.0514 


0.0553 


8.79 


4.09 


2.92 


1141.56 


531.00 


379.52 


2 


16 


0.38 


0.40 


18.5- 


-19.0 


0.0514- 


-0.0459 


0.0486 


6.50 


3.69 


2.50 


1181.82 


671.38 


454.55 


2 


12 


0.45 


0.49 



Table 4. Same as Table|3]but for the sample identified with method #2. 



Mag range 


Mass 


range 


Mid-mass 


dN 


errH 


errL 


dN/dM 


errH 


errL 


dN/dlogM 


errH 


errL 


12.0-12.5 


0.7380- 


-0.6420 


0.6900 


4.00 


3.18 


1.94 


41.67 


33.12 


20.17 


1.82 


0.58 


0.66 


12.5-13.0 


0.6420- 


-0.5750 


0.6085 


28.00 


6.36 


5.27 


417.91 


94.95 


78.62 


2.77 


0.20 


0.21 


13.0-13.5 


0.5750- 


-0.5070 


0.5410 


61.00 


8.86 


7.79 


897.06 


130.27 


114.62 


3.05 


0.14 


0.14 


13.5-14.0 


0.5070- 


-0.4200 


0.4635 


78.00 


9.87 


8.82 


896.55 


113.50 


101.35 


2.98 


0.12 


0.12 


14.0-14.5 


0.4200- 


-0.3260 


0.3730 


79.00 


9.93 


8.87 


840.43 


105.64 


94.41 


2.86 


0.12 


0.12 


14.5-15.0 


0.3260- 


-0.2440 


0.2850 


104.00 


11.23 


10.19 


1268.29 


137.01 


124.22 


2.92 


0.10 


0.10 


15.0-15.5 


0.2440- 


-0.1830 


0.2135 


97.00 


10.89 


9.84 


1590.16 


178.47 


161.25 


2.89 


0.11 


0.11 


15.5-16.0 


0.1830- 


-0.1390 


0.1610 


68.00 


9.29 


8.23 


1545.45 


211.17 


187.07 


2.76 


0.13 


0.13 


16.0-16.5 


0.1390- 


-0.1085 


0.1237 


52.00 


8.26 


7.19 


1704.92 


270.92 


235.86 


2.68 


0.15 


0.15 


16.5-17.0 


0.1085- 


-0.0869 


0.0977 


42.00 


7.54 


6.46 


1944.44 


349.00 


299.14 


2.64 


0.17 


0.17 


17.0-17.5 


0.0869- 


-0.0703 


0.0786 


26.00 


6.17 


5.07 


1566.27 


371.81 


305.69 


2.45 


0.21 


0.22 


17.5-18.0 


0.0703- 


-0.0591 


0.0647 


9.00 


4.12 


2.96 


803.57 


368.08 


264.11 


2.08 


0.38 


0.40 


18.0-18.5 


0.0591- 


-0.0514 


0.0553 


7.00 


3.78 


2.60 


909.09 


491.41 


337.41 


2.06 


0.43 


0.46 


18.5-19.0 


0.0514- 


-0.0459 


0.0486 


8.00 


3.96 


2.78 


1454.55 


719.64 


506.16 


2.21 


0.40 


0.43 


19.0-19.5 


0.0459- 


-0.0408 


0.0433 


7.00 


3.78 


2.60 


1372.55 


741.94 


509.43 


2.14 


0.43 


0.46 


19.5-20.0 


0.0408- 


-0.0369 


0.0389 


11.00 


4.43 


3.28 


2820.51 


1135.34 


840.70 


2.40 


0.34 


0.35 


20.0-20.5 


0.0369- 


-0.0331 


0.0350 


3.00 


2.94 


1.66 


789.47 


772.76 


436.40 


1.80 


0.68 


0.80 


20.5-21.0 


0.0331- 


-0.0296 


0.0314 


1.00 


2.32 


0.87 


285.71 


663.68 


247.44 


1.31 


1.20 


2.01 



the careful comparison of model isochrone s to the Pleiades 
photometric sequence by iBell et al] (|2012l ) . Other clusters 
with ages derived by the lithium depletion boundary tend 
to agree with older age estimates although with a possible 
trend towards slightly older ages, e.g. IC4665 (36 Myr vs 
28±4 Myr Mcrmilliod 1981; Manzj et al 2008) IC 2602 (30- 
67 M yr vs 46±6 Myr: IKharchenko et al.ir2005l : lOobbie et al 



20101) NGC 2547 (20-35 Myr vs 34-36 Myr: iNavlor et al 



20021 : I Jeffries fc Oliveiral l2005l) or M35 (200li^^ Myr vs 



175 M yr; ISung fc Bcsselll Il999l : iBarrado V Navascues et al.l 
l200lh . We will employ the isochrones for the lithium test 
age of 90 Myr in the case of a Per, bearing in mind the 
current uncertainty on its age of the order of 10 Myr. 

Several distance esti mates have been published for a 
Pe r: l<d0.btll Pc bv.Robichon et al.1 (1999 ) 176 2 ±5.0 pc 



bv lPinsonneault et all j 19981 ) and lMakarovl (|2006h . The lat- 
est val ue derived from a re vised reduction of the Hipparcos 
data bv lvan LeeuwenI l|2009l ) suggests a distance of 172.4±2.7 
pc which we adopt in this work. 



To s ummarise, we ado pt in this work a distance of 
172.4 pc l|van Leeuw cn 2009) and employed the Lyon group 
Next Gen l|Baraffe et al..,1998i ) and DUSTY (|Chabrier et al.1 
I2OOOI) models at an age of 90 Myr to convert the luminosity 
function into a mass function. We should point out that the 
lowest mass brown dwarfs in a Per are warmer than 1400 K, 
the upper limit wh ere the COND models should be used 
iBaraiffe et al.ll2002l) . 



7.2 The luminosity function 

In this section, we construct two luminosity functions: i) 
we used the sample of 10,176 stars in a Per with computed 
membership probabilities (Section l4.ip : and ii), the 685 can- 
didates identified with method #2 (Section 14. 2[) . The lumi- 
nosity function of the former method is derived by summing 
membership probabilities of all stars fitted to distribution 
functions in the vector point diagram, whereas the luminos- 
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Figure 9. Luminosity (left) and system mass (right) functions derived from our analysis of the UKIDSS GCS DR9 sample of member 
candidates in a Per. Error bars are Gehrels errors. The brightest bin and the last bins are very likely contaminated because of saturation 
and incompleteness, respectively. The left-hand side panel compares the luminosity function obtained from the probabilistic approach 
(black symbols and black line) and the luminosity function derived from the selection outlined by method #2 (red colour). Note that 
the sample of method #2 extends two magnitude bins fainter but they are incomplete as is the brightest bin due to saturation. The 
right-hand side panel compares the a Per mass function derived from this probabilitic approach (filled black dots linked by a solid 
line) and the mass function derived from method #2 (red symbols and red line). Error bars on the mass (x-axis) are 3(t uncertainties 
considering the errors on the age and distance of a Per. The Pleiade s mass function derived in a similar manner is overplotted in green 
for comparison along writh the field (system) mass functions in blue l|ChabrieJl2005l l. 



ity function of the latter is derived simply by summing the 
number of member candidates. 

Both luminosity functions i.e. the number of stars and 
brown dwarfs as a function of magnitude plotted per 0.5 mag 
bin is displayed in Fig. [§1 Both luminosity functions look 
very similar and match each other within the error bars. The 
numbers of objects per 0.5 mag bin increase quickly to reach 
a peak around Z — 14.5-15 and drop off afterwards down to 
the completeness of our survey with a possible peak around 
Z = 19.5-20 mag ( Tables O&ll}. The brightest bin is a lower 
limit due to the saturation limit of the GCS survey. The last 
four bins included in method #2 are not present in the prob- 
abilistic approach because the broad cluster distribution and 
low separation from the field causes the probabilities to be 
washed out. All bins in the probabilistic luminosity function 
are complete while the last two bins from method #2 are 
incomplete due to the constraints imposed on the Z-hand 
detection. Moreover, the a Per luminosity function is very 
similar to the Pleiades one deri ved in a similar man ner using 
the same homogeneous survey (|Lodieu et al.ll2012h although 
less populated mainly because of the smaller areal coverage. 



7.3 The mass function 

In this section we adopt the logarithmic f orm of th e Initia l 
Mass Function as originally proposed by ISalpeteil l|l955h : 
f(logj^Qm) = dn/dlogi(,(m) « m~^ where the exponent of 
the mass spectrum a = x + 1 following the formulation of 



I Chabrie J |200^). The Z = 12-21 mag range translates into 
masses between ~0.74 and ~ 0.03 Mq (19 mag and 0.046 
M0 in the case of the prob abilistic approach) , assuming a 
revised distance of 172.4 pc (|van Leeuwenllioog ) and an age 
of 90 Myr for which the models are computed. 

We included in Fig. [5] errors in both the x-axis (logM) 
and y-axis (dN/dlogM) as follows. For the error bars on 
the masses, we consider ed three times the uncertainties 
on the age f90 ±10 Myr: IStauffer et al.lll999l ) and distance 
(172.4±2.7 pc; Ivan LeeuwenI 120091 7 of a Per given us a va- 
lidity range of 3a on the x-axis. Hence, we computed the 
masses with the 60 Myr NextGen and DUSTY isochrones 
shifted at a distance of 164.3 pc to define the lower limit 
and repeated the procedure with the 120 Myr isochrones for 
a distance of 180.5 pc as upper limits. The uncertainties on 
the y-axis i.e. the dN/dlogM values are simply Gehrels er- 
ror bars. This a Per mass function, directly co mpared to 
the P leiades l|Lodieu et al.l [201 2) and the field l|Chabrierl 
l2005h mass functions, agree within the error bars. We should 
point out the recent mass function of the field published by 
iKroupa et al] ((2011) and described as a power-law is almost 
identical to the log-normal form of IChabrierl (|2005l ). 

In Fig. [TO] we show a log-normal fit for a Per incorpo- 
rating higher mass data points from .Prosserl (|l992l ) in order 
to provide constraint on the parameters of the fit, in partic- 
ular the characteristic mass which requires sufficient points 
on both sides of the peak in the function. We translated 
the 'corrected' luminosity function values making a small 
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AlphaPer: = 0.344±0.045; a = 0.458±0.019; = 2.275 




logi(,(m/Mg) 

Figure 10. Log-normal fit to the GCS DR9 data (metliod #2; 
triangles with error bars) in conjunction with higher mass data 
points (stars with error bars) taken from Prosser (1992). The 
least-squares fit to the data points is the solid line with the shaded 
region corresponding to a formal la uncertainty. 

update to the absolute magnitudes for the distance modu- 
lus used here (6.18) over the value of 6.1 in IProssej l|l992[ ). 
The visual band m ass-luminosity relation used comes from 
iMarigo et al.l (|2008l ) evolutionary modelfl We include in the 
fit only those high er mass points that are complete, i.e. for 
Mv < +5 from lProsser, (,1992 ), and excluding our own high- 
est mass point from the GCS luminosity function, but we 
include the four lowest mass points from the GCS since ex- 
cluding them does not significantly alter the fit. The mass 
function appears to be well represented by a log-normal with 
goodness-of-fit ~ 2.3 which indicates some systematic 
fluctuations over and above the assumed sampling errors 
that could easily be due to sample contamination and/or 
systematic errors resulting from the assumed models. 

It is interesting to compare this mass function with 
those from the Pleiades (|Lodieu et al.|[20l3 ) and Praesepe 
(Boudreault et al. subm. to MNRAS), with similar higher 
mass constraints from optical photographic plate surveys 
- see Fig. 1111 In the case of the Pleiades and Alpha Per, 
the higher mass luminosity functions have been taken as 
complete and no normalisation has been performed rela- 
tive to the GCS luminosity functions, whereas for Prae- 
sepe we found that the mass function resulting from the 
IJones fc S tauffcr ( 1991) luminosity function is discontinuous 
with the GCS mass function from Boudreault at al. (2012, 
subm to MNRAS). We determined a relative normalisation 
of 0.447 in the log (a factor 2.8) for a minimal chi-squared 
in the log-normal fit for Praesepe. 

In Table Owe compare the log-normal fit para meters to 
the fie ld sy stem mas s function parameterised by IChabrieil 
l|2003l ) and IChabrieij l|2005l ). There is some marginal evi- 
dence here for a variation in characteristic mass at the 1 
to 2(T level between a Per and Praesepe and the Pleiades, 
but this must be treated with caution given the range of 
goodness-of-fits obtained (1.0 < < 4.4) and particu- 

^ http:/ /stev. oapd.inaf.it/cgi-bin/cmd_2. 3 




-1.5 -1 -0.5 0.5 1 



loglo("»/M©) 

Figure 11. Log- normal fit to the GCS DR9 data (triangles with 
error bars) in conjunction with higher mass data points (stars 
with error bars) taken for the Pleiades (Lodieu et al. 2012a ex- 
cluding the 3 lowest mass bins; higher mass points from the 
unpublis hed compilation s of Prosser and Stauffer, see for exam- 
ple jHamblx^er^l] |l99^ and references therein); Alpha Per (this 
work); and Praese pe (Boudreault et al. su bm to MNRAS; higher 
mass points from Ijones &: Staufferl l ll99 j) ). In each case, least- 
squares fits to the data points are the solid line with the shaded 
region corresponding to a formal Icr uncertainty. 



larly the significant departure from the fit for the Pleiades 
at the low mass end. There is a clear statistically significant 
difference between the dispersion values of the field and a 
Per mass functions, not unexpected due to the difference in 
age. While we caution that the fitted values can be sensitive 
to the relative normalisation between the GCS and higher 
mass data, changes in the relative offsets tend to narrow the 
log-normal fit rather than broaden it. In any case, it is in- 
teresting to note the general log-normal trend in these wide 
mass range mass functions. 

Assuming that the observed lithiu m depletion bound- 
ary is at M ~ 0.075 Mr {Mz = 11.155; IStauffer et al.lll999l : 
iBarrado v Navascues et al]|2004h and a distance of 172.4 pc, 
the sample extracted by method #2 contains 685 a Per 
member candidates, divided up into 632 stars (92.3%) and 
53 brown dwarfs (7.7%). Lower percentages of brown dwarfs 
are obtained considering the sample of 431-728 high prob- 
ability members (p^40-60%) identified in the probabilistic 
approach, because of larger uncertainties on the probabil- 
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Population 



Characteristic 
mass mc (M©) 



Dispersion 



Alpha Per 
Pleiades 
Praesepe 



0.344 ± 0.045 
0.247 ±0.047 
0.328 ± 0.035 



0.458 ± 0.019 2.275 
0.456 ± 0.023 4.382 
0.434 ±0.015 0.962 



Field (Chabrier 2003) 0.22 
Field (Chabrier 2005) 0.25 



0.57 
0.55 



Table 5. Comparison between log-normal mass function param- 
eters for the Alpha Per, Pleiades and Praesepe clusters as de- 
termined from GCS DR9 data in conjunction with higher mass 
bin data from optical photographic proper motion surveys, com- 
pared with the field system mass fun ction parameters quoted by 
IChabriei] ||2003|') and lChabrieJ |200l). 



ities at the faint end. Hence, the star (~0.6-0.08 Mq) to 
brown dwarf (0.08-0.04 Mq) ratio in a Per spans 11.9 (10.4- 
12.7; 3(T limits using the lower and upper distance estimates) 



to 16.8 



+2.5 



_oo 0+4.6 

OO.O g , 



in agreemen t wit h measurements 



m 



IC348 (8.3-11.6; iLuhman et al] l2003l : lAndersen et al.l 

20061) but higher than other op en clusters like M35 (4.5; 
Barradq v Navascues et al. l200l | ) or the Pleia des (3.7 and 



5.7-8.8: iBouvier et al.lll998l : lLodieu et al.ll2012l '). young star- 



forming regions (3.0-6.4 for the Trapezium Cluster; 3.8-4.3 



for a Orionis; 3.8 for Chamaeleon: iHillenbrand fc Carpenter 



2000| : iMuench et a,l.l l2002l : lAndersen et all |2006|: 



Luhman 



200?!: iLodieu et al.ll2009l)T the field ( 1.7-5.3; iKroupal I2OO2I : 



Chabrier! 20051 : Andersen et al.l [20061 ). and hydrodynamical 
simulations of star clusters (3.8-5.0; ,Batc 2009, 20l2). We 
fist the ranges of the ratios because the stellar and substellar 
intervals differ slightly from study to study. 



8 SUMMARY 

We have presented the outcome of a wide (~56 square de- 
grees) and deep (J ~ 19.1 mag) survey in the a Per open 
cluster as part of the UKIDSS Galactic Clusters Survey Data 
Release 9. The main results of our study can be summarised 
as follows: 

• we recovered member candidates in a Per previously 
published and updated their membership assignations 

• we selected photometrically and astrometrically poten- 
tial a Per member candidates using two independent but 
complementary methods: the probabilistic analysis and a 
more standard method combining photometry and proper 
motion cuts 

• we investigated the if-band variability of a Per cluster 
members and found virtually no variability at the level of 
0.06-0.09 mag 

• we derived the luminosity function from both selection 
methods and found no difference within the error bars 

• we derived the a Per mass function over the 0.5-0.04 
M0 mass range: its shape is similar to the Pleiades mass 
function and best represented by a log-normal form with a 
characteristic mass of 0.34 M0 and a dispersion of 0.46. 



This paper represents a significant improvement in our 
census of the a Per low-mass and substellar population 
as well as our knowledge of the mass function across the 
hydrogen-burning limit over the entire cluster. We believe 
that this paper will represent a reference for many more 
years to come in a Per. We will now extend this study to 
other regions surveyed by the GCS to address the question of 
the universality of the mass function using an homogeneous 
set of photometric and astrometric data. Future work to con- 
strain current models of star formation includes the search 
for companions to investigate their multiplicity properties, 
the determination of the radial velocities of a Per members, 
and deeper surveys to test the theory of the fragmentation 
limit. 
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Table Al. Sample of known member candidates previously published in a Per and recovered in GCS DR9. We list the equatorial 
coordinates (J2000), GCS ZYJHK1K2 photometry, proper motions (in mas/yr) and their errors, reduced chi-squaxed statistic of the 
astrometric fit for each source (x^ value), membership probabilities when available from our probabilistic study, and names from the 

literature. A line in the probability column means that the object la^ks measurement a Per member candidates are ordered by 

increasing right ascension. This table is available electronically in the online version of the journal. 



R.A. 


Dec. 


Z 


Y 


J 


H 


K\ 


K2 


/ic«cos5±err 


± err 




Prob 


Name 


02 58 17.66 

03 01 21.38 


+48 28 00.4 
+48 35 23.3 


16.152 
13.971 


15.700 
13.664 


15.071 
13.142 


14.531 
12.494 


14.187 
12.257 


14.175 
12.267 


23.07+2.91 
24.39+2.86 


-14.86+2.91 
-21.71+2.86 


0.59 
0.11 


0.77 


DH12.Prob73.7 
DH15.Prob70.8 


03 50 37.08 
03 50 37.08 


+48 12 31.4 
+48 12 31.4 


14.124 
14.124 


13.621 
13.621 


13.079 
13.079 


12.534 
12.534 


12.234 
12.234 


12.226 
12.226 


21.45+2.03 
21.45+2.03 


-35.54+2.03 
-35.54+2.03 


6.38 
6.38 




AP265.M9.9.Y? 
DH302_Prob79.1 



Table Bl. Coordinates (J2000), near-infrared {ZY J H K1K2) photometry, and proper motions (in mas/yr) for all high probability 
(p^40%) members in a Per identified in the UKIDSS GCS DR9 using the probabilistic approach. The last column gives the membership 
probability. Sources are ordered by increasing right ascension. This table is available electronically in the online version of the journal. 



R.A. 


Dec. 


Z 


Y 


J 


H 


Kl 


K2 


IJ-aCOsS 




Prob 


02 58 52.52 
02 58 57.10 


+49 40 32.6 
+50 44 41.4 


14.543 
15.074 


14.759 


13.655 
14.213 


12.993 
13.590 


12.761 

13.335 


12.748 
13.344 


26.74 

23.90 


-22.66 
-20.86 


0.71 
0.61 


03 50 01.17 
03 50 20.08 


+48 20 57.3 
+48 13 54.8 


16.587 
15.402 


16.104 
15.029 


15.490 
14.504 


14.812 
13.940 


14.494 
13.645 


14.462 
13.617 


20.58 
25.02 


-26.20 
-19.56 


0.46 
0.43 



Table CI. Coordinates (J2000), near-infrared {ZYJHKIK2) photometry, and proper motions (in mas/yr) for all member candidates 

in a Per identified in the UKIDSS GCS DR9 with the standard method (method #2), including known members from earlier studies. 
Sources are ordered by increasing right ascension. This table is available electronically in the online version of the journal. 



R.A. 


Dec. 


Z 


Y 


J 


H 


Kl 


K2 




t^5 


02 57 51.18 
02 57 52.10 


+48 08 29.0 
+48 23 58.8 


16.810 
17.175 


16.101 
16.459 


15.536 
15.810 


14.900 
15.192 


14.598 
14.851 


14.632 
14.828 


17.56+2.96 
22.12+3.05 


-29.20+2.96 
-17.13+3.05 


03 50 18.91 
03 50 35.47 


+48 24 59.1 
+47 25 56.3 


18.766 
17.188 


17.618 
16.424 


16.684 
15.716 


16.108 
15.156 


15.610 
14.728 


15.542 
14.730 


18.12+2.54 
22.78+2.59 


-25.55+2.54 
-29.06+2.59 
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Table Dl. Coordinates (J2000), near-infrared {ZYJHKIK2) photometry, and proper motions (in mas/yr) for Y J H K-on\y (top) and 
JHK-only (bottom) detections. 





R.A. 




Dec. 




z 


Y 


J 


H 


K\ 


K2 


fXoiCOsS 




Comments 


03 


16 


26. 


.24 


+49 


nn 
00 


12. 


.2 


— 


20.314 


18.797 


17.849 


16.993 


17.030 


14.76±4.60 


-40.10±4.60 




03 


21 


14 


.97 


+49 


14 


23. 


.2 


— 


19.548 


18.220 


17.410 


16.673 


16.713 


21.93±3.59 


-18.70±3.59 




03 


23 


09 


.75 


+50 


20 


03 


.3 


— 


19.319 


18.017 


17.307 


16.591 


16.658 


18.62±4.84 


-36.35±4.84 




03 


24 


01 


.62 


+46 


48 


52 


.7 


— 


19.291 


18.038 


17.273 


16.596 


16.720 


11.99±5.96 


-14.60±5.96 




03 


z 1 


4y, 


.28 


"t-OU 


A 


zo. 


.3 


— 


18.483 


17.404 


16.771 


16.191 


16.127 


13.76±4.35 


-21.26±4.35 


detected in Z 


03 


lo 


11. 


.64 


+51 


An 


50. 


.6 


— 


18.138 


15.299 


14.980 


14.781 


14.804 


28.61±3.04 


-22.53±3.04 


detected in Z 


03 


28 


38. 


.15 


+48 


59 


51. 


.1 


— 


20.508 


18.738 


17.845 


16.997 


16.853 


26.30±4.44 


-29.45±4.44 




03 


OQ 


4y. 


.62 


1 /t 
-f-^o 


00 


Uo. 


.3 


— 


20.112 


18.739 


17.846 


16.998 


16.978 


17.76±5.21 


-35.06±5.21 




03 


QA 

oU 


oZ. 


.69 


+OU 


OQ 
ZO 


tA 
.34. 


.7 


— 


19.908 


18.498 


17.390 


16.481 


16.424 


28.15±4.34 


-30.99±4.34 




03 


QO 

oz 


z7. 


.13 


+48 


nn 
OU 


C^A 

54. 


.3 


— 


19.428 


18.138 


17.338 


16.629 


16.546 


31.10±4.60 


-23.63ib4.60 




03 


oZ 


A 


.65 


+0U 


Ui 


oy 


.8 


— 


20.449 


19.026 


17.926 


16.966 


17.130 


20.66±6.40 


-34.40±6.40 




03 


36 


03. 


.86 


+50 


39 


57. 


.7 


— 


20.269 


18.899 


17.680 


17.017 


17.059 


5.34±6.87 


-6.19ib6.87 


detected in Z 


03 


QO 

oy 


00. 


.40 


L An 
+4y 


na 
UO 


oy. 


.5 


— 


20.228 


18.785 


17.865 


17.098 


16.986 


8.58±7.88 


-22.64ib7.88 


detected in Z 


no 


59 


48, 


,00 


+47 


50 


31, 


Q 
.0 






io.olU 


i ( .\joo 


17 1 Q« 
i t . ioo 


"17 /I 

1 1 , Z4y 


A QQ-\-(i C^Q 

4.ooitD.oo 


1 A riR-UR f^Q 
— I'^l.UDztD.Oo 


no Zj,Y images 


Uo 


01 


14, 




+49 


03 


05, 


.0 






io. i\fo 


i / .yiD 


iO.OZ'^i: 


T 7 T 1 
1 ( , 1 iZ 


7 Q/1_l_fi QQ 
i .o4itD. 00 


1 n c;/i -UR QQ 
— lU.04itD. 00 


no Y image, detected in Z? 


Uo 


09 


07, 


.00 


+49 


37 


36, 


Q 
.0 






iy .UD i 


iO.Z4i 


1 7 QQO 

i / .ooy 


1 i .000 


— U.Uoity.UD 


1 1 0-Uo nfi 
— i. iyzty.Uu 


no Y image 


03 


10 


32. 


.62 


+49 


25 


19. 


.4 






19.043 


18.290 


17.289 


17.416 


U.oo±o.oo 


— o.l7±o.oo 


detected in Y 


Uo 


11 


26. 


'7fi 


+49 


13 


52. 



.z 






iy.UlO 


io.oio 


1 7 OOK 

i / .zyo 


I'? AT A 
1 / .4i4 


Q Q1 _l_0 Q1 

o.oiito.oi 


1 OR-LQ Q1 

— l.ZDiLo.oi 


uetectea m Zj+ i 


Uo 


12 


05. 


.Oi 


+49 


02 


16. 


.u 






10. y ( 


i ( .oDo 


1 7 OA A 
i / .Z44 


1 7 OQQ 
1 / .Zoo 


a 1 7-1- Q AQ 


1 QQ-i-Q A Q 
— lZ.yyit0.4y 


no Y image 


Uo 


12 


25. 


. / D 


+49 


43 


42. 


.0 






iy.Uo!$ 


io.io/ 


1 ft CiQA 

iD.yo4 


1 7 rxdft 


1 A QQ-l-7 1 
i4.00lt / .iZ 


1 R ,17-1-7 1 
— ID. 4/ it / .iZ 


no Y image 


Uo 


14 


56. 


AO 


+50 


08 


28. 


Q 
.0 






1 nA Q 


1 7 QQQ 

i ( .yoo 


17 1 00 

i r .iyz 


1 7 rifii 

1 / .UDl 


1 1 K-l-Q 01 
i. iDlto.Zi 


01 flQ-l-Q 01 

— zi.uyito.zi 


uetectea in Zj-f- 1 : ; 


Uo 


16 


22. 


.zo 


+52 


32 


00. 


Q 

.y 






io.oy i 


1 7 OQQ 

i ( .zyy 


1 K c;77 
io.O ( f 


1 R 7/1Q 
ID. ( 40 


OQ 7Q-I-7 on 

Zo. ( yit ( .zu 


1 n QQ-i-7 on 
— lU.ooit/ .zu 


spike of a bright star 


Uo 


16 


25, 


no 


+52 


32 


09, 


. i 






io.0f50 


1 7 QR1 
i / .ODi 


1 7 ni Q 
i / .Uio 


1 7 QHO 


QO t: /I _L C QC 

OZ. 04±0.oD 


OR OQ-LC^ QR 

— ZD.ZoitO.oD 


detected in 


Uo 


17 


37, 


.ol 


+47 


05 


14, 


.0 






1 Q OQQ 


1 Q 1 Qfi 

io. ioo 


17 1 A7 

i / . iu ( 


1 7 1 


1 R CrO-LQ /I Q 

lD.0y±o.4o 


7 QQ_I-Q A Q 

— / .vSoito.4o 


no Y image 


Uo 


17 


49, 




+46 


58 


35, 


Q 
.0 






1 o nc;Q 
iy.U0f5 


1 Q 1 
io. ioD 


1 7 onn 
i / .zyu 


1 7 oon 
1 i .ZzU 


AOJ_n OQ 
O.ozity.zo 


Q n7-l-0 OQ 

— o.y / zty.zo 


detected in Y 


Uo 


18 


23. 


.yt) 


+46 


26 


49, 


.0 






io.O / 


1 7 f^t^Q 

i / .OOo 


iu.oyy 


1 R 77/1 

Id. m 4 


D.oydiO.Oi 


AO Q/1-LK c;i 

— 4z.o4itO.Oi 


detected in Y 


Uo 


19 


11. 


no 


+51 


24 


47. 


.U 






ly.uyo 


1 /IRQ 
i0.4DO 


1 7 Q1 A 

i / .oi4 


1 7 Qon 

1 i .ozy 


oo-Lo c;i 
z.Zzity.Oi 


— y.ou±y.oi 


no Y image 


Uo 


19 


19. 


1 o 

.ly 


+46 


10 


18. 


Q 
.0 






ly.uui 


1 7 7^;q 


1 ^; noi 
io.yzi 


iD.y / 


OQ ccn-Ln CR 
zo.oUiLy.OD 


1 Q CO-Un KR 

— io.ozity.oD 


aetectea m Z/+ 1 


Uo 


20 


41. 


. / y 


+50 


45 


38. 


f! 
.0 






io. / yo 


1 43 1 Q7 

io.io/ 


1 7 riQn 
i ( .UoU 


1 7 /I Q/t 
i / .4o4 


— U.yDity.4D 


1 n ACi-i-Ci Ad 

— iu.4yity.4D 




Uo 


21 


14. 


'7A 
. 1 4 


+46 


36 


27. 


1 
.± 






iy.UU/ 


1 Q 1 ,1 Q 
iO.i40 


17 1 on 
i r .ioU 


1 7 oon 
i / .zzU 


Q nKJ-O QQ 

o.yoity.oo 


1 Q KQJ-O QQ 

— io.ooity.oo 


uetectea in Zj+ i i : 


Uo 


21 


53. 


AA 


+46 


47 


02. 


.D 






iy.uyz 


iO.ZUi 


1 7 QC;7 
i ( .00 ( 


17 1 QQ 


a Al -\-m AQ 


1 QQ-l-1 n AQ 
— iZ.oyitlU.4o 




Uo 


23 


02, 


1 /I 


+52 


13 


58, 


Q 
.0 






io.y4o 


1 7 Q7Q 
i / . / 


1 7 nnn 
i / .UUU 


17 1 on 
1 i .IzU 


1 1 7nj_ in c; /I 
li. /U±iU.04 


R 7R 1 1 n K /I 
— D. (DltiU.04 


detected in Y 


Uo 


24 


03, 


.U / 


+50 


03 


01, 


.U 






1 Q HOC 

iy.uyo 


1 7 OT 1 

i / .yii 


1 7 n7i 
i / .U / i 


1 7 non 
1 i .uyu 


1 ft Q c;_l_7 OQ 

ID.oOit i .yo 


on c;0-l-7 OQ 

— zu.ozit ( .yo 


detected in Y 


Uo 


24 


32, 


.UU 


+47 


04 


29, 


.0 






1 Q Rc;r 
iO.DOD 


1 7 O/l 1 

i / .y4i 


1 7 HQ 
i / .UoO 


1 R Q7Q 

ID.o i 


OO-LQ no 
O.zyito.Uz 


1 1 AO-LQ no 
— li.4zito.Uz 


uetectea m /*+ 1 


Uo 


24 


46, 


OA 

.z4 


+46 


36 


25, 


A 
.4 






iy.U4U 


1 Q OQO 

io.zoy 


1 7 O/l /I 
i / .z44 


1 7 1 

1 i .IDO 


1 n 7c;_L 1 n Qn 
lU. /OitiU.oy 


01 OO-LI n QO 

— zi.zZitiu.oy 




03 


26 


14, 


.34 


+51 


55 


36, 


.4 








18.718 


17.778 


16.824 


16.717 


2.88±8.33 


-7.91+8.33 




03 


27 


14. 


.93 


+52 


15 


58. 


.3 






18.750 


17.710 


16.717 


16.713 


6.39ib6.88 


-9.24+6.88 


detected in Z+Y 


03 


27 


32. 


.27 


+47 


11 


45. 


.4 






19.079 


18.325 


17.133 


17.227 


-1.16ibl0.13 


-16.30+10.13 




03 


27 


43. 


.73 


+46 


55 


02. 


.9 






19.042 


18.433 


17.275 


17.368 


5.02itl0.04 


-8.77+10.04 




03 


28 


16. 


.47 


+48 


29 


41. 


.9 






18.986 


18.067 


17.075 


16.986 


9.54ib6.66 


-10.33+6.66 


detected in Z+Y 


03 


30 


17, 


.49 


+48 


04 


56, 


.8 






19.040 


18.184 


17.283 


17.416 


13.82±7.76 


-30.29+7.76 




03 


30 


46, 


.17 


+45 


57 


36, 


.0 






18.938 


17.983 


17.221 


17.139 


9.57±8.37 


-9.08+8.37 


detected in Z+Y 


03 


30 


55, 


.78 


+45 


55 


56, 


.6 






18.465 


17.740 


16.968 


16.958 


25.29±8.24 


-19.76+8.24 


detected in Z+Y 


03 


31 


01, 


.32 


+46 


09 


14, 


.4 






19.092 


17.954 


17.266 


17.256 


2.25±9.52 


-0.70+9.52 


detected in Z+Y 


03 


31 


08, 


.17 


+50 


10 


16, 


.6 






18.924 


18.075 


17.172 


17.177 


-2.82±9.25 


-5.09+9.25 




03 


34 


53. 


.62 


+47 


34 


24. 


.5 






19.077 


18.227 


17.341 


17.419 


20.62ib9.37 


-4.52+9.37 


detected in Z+Y 


03 


36 


11. 


.61 


+46 


48 


35. 


.0 






18.925 


17.920 


16.965 


16.953 


14.27ib7.87 


-21.11+7.87 


detected in Z+Y 


03 


36 


26. 


.40 


+48 


38 


22. 


.4 






19.081 


18.295 


17.183 


17.317 


-2.90ib9.79 


0.75+9.79 


detected in Y 


03 


43 


15. 


.74 


+47 


34 


45. 


.0 






18.859 


17.886 


17.177 


17.267 


22.76ib8.90 


-10.89+8.90 


detected in Z+Y 



18 N. Lodieu et al. 



Table El. Coordinates (J2000), near-infrared {ZYJHKIK2) photometry, and proper motions (in mas/yr) for substellar multiple system 
candidates identified photmetrically in a Per 
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